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A tetragonal (Bain path) distortion of a metal with an fee (bee) ground state will initially cause 
an increase in energy, but at some point along the Bain path the energy will again decrease until a 
local minimum is reached. Using a combination of parametrized tight-binding and first-principles 
LAPW calculations we show that this local minimum is unstable with respect to an elastic distortion, 
except in the rare case that the minimum is at the bec (fee) point on the Bain path. This shows 
that body-centered tetragonal phases of these materials, which have been seen in epitaxially grown 
thin films, must be stabilized by the substrate and cannot be free-standing films. 



PACS numbers: 61.66.-f,64.60.-i,68.55.-a 

I. INTRODUCTION 

The concept of the tetragonal distortion of a face- 
centered cubic lattice into a body-centered cubic lattice, 
as well as the reverse transformation, was introduced by 
Bain in 1924i as a model for the formation of Martensite 
in steels. Although this model is not exactly correct^ 
the concept of the "Bain path" has been useful in inves- 
tigating the behavior of element and compound metallic 
systems. This usefulness has increased with the develop- 
ment of total energy electronic structure methods, which 
can investigate regions of the Bain path not accessible by 
experimental methods. By following the Bain path from 
fee to bec, Mehl and Boyer— (in the following referred 
to as MB) showed that the fee elements aluminum and 
iridium are elastically unstable in the bec phase, i.e., the 
Bain path reaches a local maximum at the position of the 
bec lattice. Wills, et al.£ did first-principles calculations 
of the Bain path of gold and platinum which showed that 
these elements were also unstable in the bec phase. In ad- 
dition they used the reverse Bain path, from bee to fee, to 
show that tantalum and tungsten were elastically unsta- 
ble in the fee phase. Craievich, et al.^^ looked at the fee 
elements nickel and rhodium, which were unstable in the 
bec state, and the bec elements cobalt, chromium, molyb- 
denum, niobium, tantalum, and vanadium, which were 
elastically unstable in the fee state. Jona and Marcus 
showed that both bec coppei— and palladium 8 , are also 
elastically unstable along the Bain path. Even metals 
which arc only metastable in the fee phase, e.g., cobalt, 
manganese, osmium, rhenium, ruthenium, titanium, and 
hafniumj£*2*i2i are elastically unstable to Bain path dis- 
tortions in the bec phase. 

The above lists demonstrate that for elements to the 
right of the alkaline earths on the periodic table, a sta- 
ble or metastable fee (bec) phase implies an elastically 
unstable bec (fee) phase. This, however, does not end 
the search for metastable phases along the Bain path. 
As noted by Craievich, et al.^ the existence of a min- 



imum and a local maximum in the Bain path implies 
the existence of another local minimum. This minimum 
must be outside the classical Bain path from fee to bec, 
though this does not preclude the existence of yet an- 
other minimum between fee and bec. ME 3 found a lo- 
cal minimum in the extended Bain path of aluminum by 
compressing along the [001] direction until c/a = 0.567, 
where c/a = 1 is the fee phase and c/a = 1/^/2 the 
bec phase. Iridium was found to have a minimum at 
c/a = 0.566. These minima are very close to the spe- 
cial value c/a — 1/V3 ~ 0.577, where the body-centered 
tetragonal lattice is 10- fold coordinated (bctio). Simi- 
lar minima are found for copper gold^ palladium, 8 
and fee titaniumAi^ This crystal structure has been ob- 
served in the a phase of protactinium and in the phase 
of mercury^ and was given the Strukturbericht designa- 
tion A a p!£ 

For bec crystals, the local minimum is at c/a > ^/2, 
where c/a = 1 is the bec lattice and c/a = \J~2 is the fee 
lattice. Sliwko, et al.^ found a local minimum in vana- 
dium at c/a ~ 1.8. There is no increased coordination 
at this site. Experimentally this structure is observed in 
indiumjii and has the Strukturbericht designation A6i^> 

Given the increased coordination near the secondary 
minimum in the Bain path for fee crystals, one is tempted 
to conclude that this state is metastable. This turns out 
not to be the case. As found by Boyer^ there is an or- 
thorhombic "magic strain" (properly, Physically Allowed 
Lattice-Invariant (PALI) strain) which transforms one 
fee lattice into another fee lattice. A similar orthorhom- 
bic magic strain transforms one bec lattice into another. 
MB showed that the magic strain path for fee aluminum 
passes through the secondary minimum of the aluminum 
Bain path, and that, when viewed in the space of the or- 
thorhombic strains, this point is actually a saddle point 
and hence unstable. In terms of the elastic constants, 
Cn — C\2 for the secondary minimum is negative. Jona 
and Marcus demonstrated this instability by direct calcu- 
lation of the elastic constants of copper 7 and palladium 8 
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at the secondary Bain path minimum. 

The concept of magic strains was generalized by Van 
de Waalii and applied to the study of metastable struc- 
tures in various systems, including silicon*^ and model 
two dimensional lattices^ In this paper we use the for- 
malism to show that the secondary minimum along the 
Bain path of an elemental fee or bcc solid is unstable 
with respect to an orthorhombic shear. We note that this 
need not be the case for more complicated unit cells - the 
metastable BCT5 state of siliconizes was found using the 
concept of magic strains. The magic strain formalism is 
a useful tool for searching for metastable structures^ 
Indeed, this paper can be considered to be a search for 
metastable states due to orthorhombic distortions of sta- 
ble fee and bcc crystals. In this case, however, we will 
show that there are no such states. 

In Section [HJ we show an extension of the Bain path 
for both fee and bcc metals. In Section ITTTl we show the 
construction of the orthorhombic strains required for fee 
— ► fee and bcc — > bcc magic strains and explore large 
regions of strain space using parametrized fits to first- 
principles results via the NRL Tight-Binding method^ 
These results show that, within the tight-binding model 
of fee and bcc metals the Bain path minima discussed 
above are elastically unstable. In Section [Tvl we employ 
first-principles density functional theory methods to di- 
rectly compute the elastic constants for structures at the 
secondary Bain path minimum. These more accurate re- 
sults show that such structures are indeed elastically un- 
stable. We conclude in Section Ivl with a brief summary 
and discussion of our results. 



II. EXTENDING THE BAIN PATH 

We can represent the Bain path of an elemental com- 
pound by the standard body-centered tetragonal unit 
cell, 
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If c = a in (JIJ we have a bcc lattice, while if c = %/2a we 
have an fee lattice. Alternatively, we can represent the 
Bain path using a /ace-centered cubic unit cell, 
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where 



a' = V2a 



(2) 



(3) 



c = a! I y/2 we have a bcc lattice. We stress that Q and 
are different representations of the same lattice. In 
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FIG. 1: Bain path for platinum as a function of c/a in 
using the tight-binding parameters described in the text. The 
volume of the unit cell was chosen to minimize the total en- 
ergy at each value of c/a. The lines indicate the positions 
of the bet io, bcc, and fee states. The top x-axis shows the 
corresponding value of c/a' from (|2J. 



this case, the unit cell represented by (JSJ) is related to the 
one represented by QJ by a rotation of 45° about the z 
axis. 

To illustrate our techniques we begin with platinum, 
which has an equilibrium fee structure. Fig. ^ shows the 
energy along the platinum Bain path as a function of 
c/a. The total energy of the system was evaluated us- 
ing the NRL tight-binding parameters for platinum22i2». 
The Brillouin zone was sampled by using a regular mesh 
of fc-points generated at c/a = 1, i.e., the bcc lattice, 
but using the body-centered tetragonal symmetry of the 
system (space group 14/mtnm — Dg^). This turns out to 
be the most efficient way to generate fc-points for all of 
the strains covered in this calculation. It does, however, 
require a large number of mesh points. We used 2119 
points in the irreducible part of the Brillouin zone, and 
carefully checked our results using other k-point meshes 



to insure convergence^* The total energy of the system 
was computed using the method described by Gillan 25 , 
with a Fermi broadening temperature of 5 mRy. By com- 
paring with smaller k-point meshes, we estimate that the 
numerical error in the calculation of the total energy was 
less than 0.1 mRy for c/a < 1.5, and less than 0.2 mRy 
over the entire range of Fig. ^ 

As expected, the global minimum in Figgis at the fee 
position, c/a = y2. (Obtaining this result was consid- 
ered a necessary condition to show k-point convergence.) 
As for other fee elements***!* there is a local minimum 
near the bctio structure, which is at c/a = \/2/3 ~ 0.816 
using representation Q). Both protactinium and the {3 
phase of mercury exist in stable or metastable configura- 
tions near this bctio structured* 

In a similar way we can study the extension of the 
Bain path for an element with a bcc ground state. This is 
most easily done starting with the face-centered primitive 
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FIG. 2: Bain path for tungsten using the tight-binding pa- 
rameters and fc-points described in the text. The notation is 
the same as in Fig. 




-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 



vectors ©, constructing a regular k-point mesh at c = a', 
(the fee lattice), again using space group I4/mmm — 
D^ h . We use a mesh with 4069 points in the irreducible 
part of the Brillouin zone. Figure |2l shows our results 
for tungsten, using our tight-binding parametersiS2i2& We 
note that there is a local minimum near c/a' = 1.15. This 
minimum is required by the continuity of E(c/a), since 
we know that the fee structure of tungsten is elastically 
unstable and so a local maximum. Unlike the case of 
the strained fee elements, there is no high-coordination 
structure near this point. Interestingly, we do see a small 
shoulder in the energy-strain curve at the location of the 
bctio structure. Experimentally, Indium exists in this 
structure, with c/a' = 1.084^ 



III. ORTHORHOMBIC STRAINS 

In the previous section we saw that the extension of 
the Bain path has a secondary minimum for elements 
with either fee or bec ground states. It is natural to 
ask whether or not this secondary minimum is in fact 
metastable. Here we show a procedure, based on MB, 3 
which answers this question, and, in the process shows 
that there is no evidence such a metastable state exists 
for any element at zero temperature. 



A. Elements with an fee ground state 

The first check in determining the metastability of 
the bet structure in Fig. ^ would be to look at long- 
wavelength distortions of the unit cell, i.e., the elastic 
constants. We will, in fact, use the technique outlined in 
MB which will allow us to look at general strains about 
the bet position. 

We begin by writing the primitive vectors for an ar- 
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FIG. 3: Map of the total energy of the lattice described by 
@ for platinum, using the NRL tight-binding parameters, 
and choosing the minimum energy volume at each point. The 
heavy straight lines represent the Bain path from bec (at the 
origin) to one of three fee lattices. The lighter straight lines 
represent the extended Bain path. The labels indicate the 
direction of the tetragonal lattice along the Bain paths. The 
fee lattice points are the minima near the labels. The bec 
lattice point, at the origin, is 10 mRy above the energy of the 
fee lattice, and the contour interval is 1 mRy. The coordinates 
(r,9) are described in J^J. 

bitrary body-centered orthorhombic unit cell in a form 
reminiscent of 




where V is the volume of the unit cell, and 

x = exp[r cos(6* + 27r/3)] 

y = exp[r cos(9 — 2-7r/3)] , and 

z = exp[rcos(9] , (5) 

so that xyz — 1 for all values of (r, 9) . We immediately 
see that the energy of this unit cell, E(V,r, 8), is peri- 
odic in 9 with period 27r/3, and that r = represents 
the bec lattice for all values of 9. Furthermore, in the 
special cases 9 = (0, 2-7r/3, — 27r/3) the orthorhombic lat- 
tice becomes tetragonal, with the tetragonal axis in the 
z, x, or y direction, respectively. These lines are just the 
extended Bain path of the system, stretched along the 
[001], [100], or [010] direction. We find the face-centered 
cubic lattice when r = (ln2)/3 along any of these paths. 
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We used our platinum tight-binding parameters to map 
the energy surface E(V,r,9) of minimizing the total 
energy with respect to volume at each point (r, 9) . The 
space group is now Immm— _D|f (^71), and we used a 
mesh of 3925 fc-points in the irreducible Brillouin zone. 
The results are shown in Fig. 01 In particular, we note 
that along each Bain path, there is a saddle point about 
5 mRy above the fee ground state. These points repre- 
sent the same lattice as the local minimum structure in 
Fig. We see that this point is unstable with respect to 
an orthogonal distortion, which means that a least one 
linear combination of the lattice constants of the tetrag- 
onal structure must be negative. To see which one, let 
ao and Co be the lattice constants representing the sec- 



ondary minimum in the Bain path using the primitive 
vectors {Q. For a crystal with an fee ground state and 
a maximum at the bec position, we will have co/ao < 1, 
though this is not necessary for the following analysis. If 
we use Q and JSJ to represent this structure we have 

V =V = 1/2 a 2 c , 
r = ro = 2/3 lnco/ao , and 
9 = 9 =0. (6) 

Let us fix the volume and expand r = ro + 2a and 9 about 
the equilibrium point. Then to first order in a and 9, Q 
becomes 



ai \ / -|a (l + 
a 3 J \ la {l-a + [3) 



ko(l 



a (l 



a 
a 



■0) 
ft 



co(l 
co(l 



2a) 
2a) 



x 
V 



(7) 
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where 

13 = V3r 9/2 . (8) 

In terms of the standard elastic strain parameters for the 
body-centered tetragonal lattice we have 

e\ = —a + (3 , e 2 = —a — j3 , e$ = 2a , 

e 4 = e 5 = e 6 = . (9) 

The lowest order contribution to the elastic energy asso- 
ciated with this strain is 

AE = V[(C n + C 12 - 4C 13 + 2C 33 )a 2 + (C u - C 12 )/3 2 } . 

(10) 

If this structure was stable, the coefficients of the a 2 and 
1 terms should be positive, in agreement with Born and 
Huang's derivation of the elastic stability criteriai^i We 
note that keeping 9 — and varying a moves the lattice 
along the Bain path. Since the point a = is a minimum 
along that path, the first linear combination of elastic 
constants must be positive. However, as seen in Fig. |3 
fixing a = and varying 9 (and therefore (3) lowers the 
energy, implying that Cu — Cn < 0. 

B. Elements with a bec ground state 

Our examination of the stability of the bet structure 
(Fig. |2J of an element with a bec ground state is similar 
to the above discussion, with minor differences because 
we are now describing a state with c/a expanded from 
the equilibrium, rather than compressed as it is for fee 
elements. 

We therefore begin with an orthorhombic distortion of 
the fcc-like primitive vectors |(2} : 



where V is the volume of the unit cell and 

x = exp[— r cos{9 + 2ir/3)] 

y = cxp[— r cos(0 — 27r/3)] , and 

z = exp[— rcos9] . (12) 

In analogy with (0} , there are three Bain paths along the 
directions 9 = (0, ±27r/3), but now r = is an fee lattice, 
and at the special value r = (In 2)/3 along any one of the 
Bain paths we find a bec lattice. 

Using the tungsten tight-binding parameters, we 
choose to describe the space group for l|12l) as 
Fmmm— D|?(^69). We use a mesh of 7813 fc-points 
in the irreducible part of the Brillouin zone. The re- 
sults of the calculation are shown in Fig. 21 The saddle 
points located along the extended Bain path, approxi- 
mately 30 mRy above the ground state, represent the 
same structure as the local minimum found in Fig. [21 We 
can determine the nature of the elastic instability here in 
a manner similar to Sec. 1111 AI Let a' and cq specify the 
local minimum on the Bain path in Figs. |2 and 01 using 
the primitive vectors @. In terms of Ijllfl . this occurs at 

V =V =l/4a 2 c , 
r = ro = — 2/3 1nco/a , and 
9 =9 Q = . (13) 

Taking r — r + 2a and expanding the lattice parameters 
in powers of a and 9, we find that (TTJ 

can be written as 
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FIG. 4: Map of the total energy of the lattice for tung- 
sten, using the NRL tight-binding parameters, choosing the 
minimum energy volume at each point. The origin of the 
graph is the fee lattice, the heavy straight lines radiating from 
the lattice are the Bain paths, and the lighter straight lines 
are the extensions of the Bain path. The labels indicate the 
orientation of the tetragonal axis along each Bain path. The 
contour interval is 2 mRy, and the energy of the fee struc- 
ture is 35.8 mRy above the energy of the bec structure. The 
coordinates (r, 8) are described in Q12|l . 

an orthorhombic distortion of J2J- To lowest order, the 
clastic strain parameters are 

ei = a + , e2 = a — (3 , e 3 = -2a , 

e 4 = e 5 = e 6 = , (14) 

where j3 is denned in (JSJ) . Thus the elastic energy associ- 
ated with this strain is identical to (|10fl , and, once again, 
we see that at the Bain path minimum Cn — C12 < 0. 

Traditionally, elastic constants in centered-tetragonal 
lattices are computed using the body-centered form of the 
primitive vectors Q rather than the face-centered form 
(J2J. As we noted, this corresponds to a 45° rotation of 
the coordinates around the z axis. Under this rotation 
the elastic constants transform so tha1s2& 
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When we compute energy change of the strain (|14(l on 
the structure defined by (|13|) using the body-centered 
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FIG. 5: Total energy as a function of c/a at the equilibrium 
volume for elements with a bec ground state. 

tetragonal primitive vectors Q), then, the energy changes 
by 

AE = V[(Cn + C 12 - 4Ci3 + 2C 33 )a 2 + 2C 66 /5 2 ] . (16) 

The saddle point behavior of the energy surface near this 
minimum leads to the prediction that Cqq < 0. 

IV. ELASTIC CONSTANTS OF THE 
STRUCTURES AT THE BAIN LOCAL MINIMA 

In the previous section we showed that for typical fee 
and bec structured elements the secondary minima in the 
Bain path, as shown in Figs.n an d|21 are actually saddle 
points in the space defined by the orthorhombic unit cells 
@J and (fTTfl . From Figs. |3and^ we see that the only 
elastically stable structures along the Bain path are the 
fee and bec lattices. In rare cases both of these struc- 
tures may be minima on the Bain path. In this section 
we will check this conjecture via first principles density 
functional theory calculations using the Linearized Aug- 
mented Plane Wave (LAPW) method^ as implemented 
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TABLE I: Position of the secondary minimum along the ex- 
tended Bain path, c/a, for the elemental metals which have 
fee or bee ground states. We also show the energy difference 
between the phases. All calculations were done by minimizing 
the total energy as a function of volume for fixed c/a, using 
the primitive vectors in Q. Energy differences are in mRy. 
Note that for strontium the secondary minimum coincides 
with the bec lattice, and for barium the minimum coincides 
with the fee lattice. 

Metal Equil. c/a \E hcc - Ef cc \ \E hcc - E bct \ \E tcc - E h ct 
Al fee 0.8f4 7.287 f.544 5.743 

Ag fee 0.894 2.404 0.f64 2.240 

Au fee 0.852 f.785 0.307 f.478 

Ba bee 1.414 f.060 L060 0.000 

Ca fee 0.911 f.648 0.f75 f.473 

Cu fee 0.935 2.706 0.033 2.673 

Ir fee 0.801 48.997 23.684 25.313 

Mo bee 1.765 30.301 23.450 5.851 

Nb bee 1.787 24.186 10.245 13.941 

Pd fee 0.876 3.752 0.405 3.347 

Rh fee 0.814 26.468 10.644 15.824 

Sr fee 0.707 0.730 0.000 0.730 

V bee 1.799 18.330 7.332 10.998 

W bee 1.687 36.471 29.195 7.276 



TABLE II: Elastic constants of at the secondary minimum in 
the Bain path for the elements described in TableQ] with C' = 
(Cu — Ci2)/2. All elastic constants are computed using the 
orientation and unit cell described by Q, and are expressed 
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proximation (GGA) exchange correlation functional of 
Perdew, Burke, and Ernzerhofi^i 

We begin by determining the position of the secondary 
minimum in the Bain path for cubic elemental metals. 
In Fig. we show the results for several of the elemental 
bec metals. In all cases there are two minima in the 
calculation, one at the bec site and one at or beyond the 
fee position on the Bain path. The results are shown 
in Table [1 where the ratio c/a is determined using the 
representation of the body-centered tetragonal lattice in 
(JIJ . With the exceptions of barium and strontium, all of 



the metals shown here have a secondary minimum of the 
type described in Sec. [H] if the ground state is fee, the 
minimum is at c/a < 1, while if the ground state is bec, 
the minimum is at c/a > s/2. Barium and strontium are 
unique in that the secondary minimum is at the fee point 
on the Bain path for barium, as shown in Fig. EI an d the 
at the bec point for strontium. We will not consider these 
elements further in this paper. 

For the remaining elements we compute the elas- 
tic constants of the tetragonal system using standard 
techniques, 7 ' 28 using the primitive vectors (QJ to describe 
the lattice. We present the results of our calculations in 
Table [H] As predicted, all of the fee ground state ele- 
ments have C = {C\\ — Cyi)/2 < at the secondary min- 
imum, and all bec ground state elements have 6*66 < 0, 
in agreement with the predictions of Sec. IIIII 

V. SUMMARY 

As noted in the introduction, it has previously been 
shown that the secondary minimum is elastically unsta- 
ble in aluminum, 3 copper, 7 iridium^ and palladium^ To 
our best knowledge the stability or instability of the sec- 
ondary minimum of a bec element has never before been 
determined. The calculations done here show that if the 
energy along the Bain path for a metal is at a minimum 
at the fee (bec) structure and a maximum at the bec (fee) 
structure, then there is a secondary minimum with a bet 
structure with is elastically unstable. 

The implication of this work is that there can be no free 
standing body-centered tetragonal structures for these 
metals. Any bet structures observed, such as in thin 
films2£ must be stabilized by an external force. In the 
case of epitaxially grown thin films, this is force is pre- 
sumably supplied by the substrate. For sufficiently thick 
films this supporting force will be insufficient to stabilize 
the bet structure, and the film will revert to its ground 
state, as Bencok et alm& observed for thin vanadium films 
grown on Fe(100). 

While we have not discussed the case of hep met- 
als such as titanium, where epitaxial bet structures are 
seenj2i24 we note that the minimum energy bec structure 
of Ti is known to be elastically unstable^ The conclu- 
sions of this paper thus apply to Ti as well. 
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